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We propose a scheme to generate and detect various kinds 
of quantum entanglement in a spin-1 Bose-Einstein conden- 
sate. It is shown that substantial many-particle entanglement 
can be generated directly in the spin-1 condensate by free dy- 
namical evolution with a properly prepared initial state. The 
scheme also provides a simple method to generate three-mode 
entanglement in the second-quantization picture and to detect 
the continuous variable type of entanglement between two ef- 
fective modes in the spin-1 condensate. 

The achievement of atomic Bose-Einstein condensate 
(BEC) [1] has raised extensive experimental and theoret- 
ical studies in this area, and these studies have opened 
up some possibilities for important applications. An im- 
pressive example of these applications is to use BECs 
with internal degrees of freedom for generation of quan- 
tum entanglement [2-7], which is the essential ingredient 
of many quantum information protocols [8]. To gener- 
ate entanglement, one requires an experimental system 
which can be prepared in a pure state, and there should 
be significant interactions between the relevant particles 
which can be coherently controlled. BECs fulfil these re- 
quirements, and therefore they could provide an ideal ex- 
perimental system for studying quantum entanglement. 
It has been shown that with coherent collision interac- 
tions in BECs, it is possible to generate substantial many- 
particle entanglement in a two-component condensate [2] , 
or to generate entangled atomic beams from the outputs 
of a single condensate prepared in a specific internal level 
which enables spin-exchange collisions [3,4]. 

A good experimental system for demonstration of these 
proposals is the Sodium condensate confined in an optical 
trap [9-11]. Recently, Stamper-Kurn et al. have realized 
an optically trapped BEC in which all the three ground 
Zeeman levels of sodium atoms are involved [9] . In such a 
three-component (so-called spin-1) condensate, there are 
spin-exchange collisions, which constantly mix the differ- 
ent spin components, and other kinds of internal spin dy- 
namics, which result in complex nonlinear collision phase 
shifts [12-17]. In this paper, we will investigate whether 
one can directly prepare and detect quantum entangle- 
ment in such a spin-1 condensate. The motivation of this 
work is two-folds: firstly, direct detection of quantum en- 
tanglement in the realized spin-1 condensate eliminates 
the requirements of the level shifting and the output tech- 
niques in previous works [2-4], and thus allows for a eas- 
ier experiment; and secondly, the spin-1 condensate with 
the presence of both the spin-exchange collisions and the 
nonlinear phase shifts is a more complicated system, and 
opens up possibilities for richer physical behaviors. We 



will show that various kinds of entanglements are avail- 
able in this system. We can get substantial many-particle 
entanglement (the spin squeezing type of entanglement 
[18,19]) as well as three-mode and two-mode entangle- 
ment (the continuous- variable type of entanglement [20]). 
All these entanglements are useful for different purposes 
of applications. 

The spin-1 BEC has been studied both theoretically 
[12-17] and experimentally [9-11], with the emphasis on 
the ground-state structure, the spin-mixing dynamics, 
and the response to the external magnetic field. Here, we 
will focus on its entanglement property. In Ref. [15], it 
was shown as a by-product that the spin ground state of 
the sodium condensate with the anti-ferromagnetic inter- 
action should be entangled. Unfortunately, this ground- 
state entanglement was not confirmed by the experiment 
[10] because it is very hard to cool the system to the spin 
ground state (the energy gap between the spin ground 
and the first excited states is many orders smaller than 
the chemical potential) , and the spin ground state is very 
fragile under the influence of the external magnetic field 
[17]. Here, we describe a different method to generate 
and detect quantum entanglement by free dynamical evo- 
lution with a proper initial state. As we will see, this 
provides a feasible method of detecting entanglement in 
the spin-1 BEC with the current technology. 

Before we enter the discussion of the spin-1 conden- 
sate, let us first clarify different kinds of entanglements 
available in a bosonic many-particle system. Consider N 
bosonic identical particles with each particle having three 
internal levels |-l-l), | — 1), and |0). First, the N particles 
are said to be entangled in the first quantization picture 
if their internal density matrix p cannot be decomposed 
into 
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where the coefficients pk are positive real numbers sat- 
isfying X^fc^fe ~ ^' ^^'^ pf is a 3 X 3 density matrix of 
the ith. particle. To check whether an Ai"-particle sys- 
tem is entangled, a good experimentally detectable crite- 
rion is the squeezing parameter. For particles with only 
two internal levels \a} and \b), one can construct a col- 
lective spin operator Jab = J2i=iJi^ with the individ- 
ual spin-1/2 operators jf = (|a)^ (6| + \b). (a|)/2, jf = 
*(|fl).(H-|fe),H)A and jf = (|6),(6|-|a),(a|)/2. 
The squeezing parameter ^"^^ is defined directly from the 
variance and mean values of the Jab operator 
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where J"{, 



n • Jab, and the ns are niutuaUy orthogonal 
unit vectors. For two-level particles, the spin squeezing 
parameter (2) with a value less than 1 provides a suffi- 
cient condition for many-particle entanglement [2] (such 
a case is also referred as spin squeezing [18,19]). For 
three-level particles, we can construct two orthonormal 
states I a) and \h) from arbitrary superpositions of the 
three levels |-l~l), | — 1), and |0), and define the lab and 
^"j, in the same way. Then, ^"^ < 1 assures that the 
projected A^-particle density matrix PpP onto the sub- 
space spanned by the states \a)^ and \b)^ (i = 1, 2, • • • , N) 
is entangled. Since the projection P is a local operator 
which cannot increase entanglement, the original density 
matrix p for the A^ three-level particles should also be en- 
tangled if ^"j, < 1. For our later applications, we define 
the spin operators J± and the squeezing parameters ^^ 
respectively by choosing \a) — |±) — (|+1) ± | — 1)) /\/2, 
\b) = |0), and ni = [cos 0, sin 6*, 0] in Eq. (2). Either of 
^^ less than 1 for some is a signature of many particle 
entanglement. Note that the spin squeezing type of en- 
tanglement represented by a small ^"f, is useful in practice 
since the inverse of ^"j, gives exactly the improvement of 
the single-to-noise ratio for atomic clocks relative to the 
standard quantum limit with the Ramsey scheme involv- 
ing the levels \a) and \b) [19]. 

The above system can also be described in the sec- 
ond quantization picture. We have three bosonic modes 
o-a {oi ~ 0, ±1) corresponding to the three internal levels. 
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modes can also be entangled like the individual parti- 
cles. Similar to the definition in Eq. (1), we say that 
two modes Cq , ap are entangled if their density matrix 
cannot be expressed as ^kPkPa ® p|j ^-nd all the three 
modes aa,ai3,a^ are intrinsically entangled if their den- 
sity matrix cannot be expressed as '^^PkPa'^ Pa-yt where 
a, /3, 7 = 0, ±1, pk are positive real numbers, and p^/^ ^ s 
is the a density matrix of the corresponding mode(s). 
Note that the mode entanglement in the second quanti- 
zation and the particle entanglement in the first quanti- 
zation are two different concepts, and they are useful for 
different purposes of applications. One can easily find 
example states which are mode entangled but not par- 
ticle entangled, or particle entangled but not mode en- 
tangled [22] . Mode entanglement of bosonic atoms could 
be useful in various continuous variable or qubit quan- 
tum information processing schemes currently based on 
entangled optical modes [8]. For two modes a±i, there is 
an easy experimentally detectable entanglement criterion 
[20]. One can define quadrature phase amplitudes X^^ = 

«a^i + e-^ Vi) /\/2 and X« = {Xl,±Xt,)/V2 
for the two modes ±1. It has been shown in Ref. [20] 
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provides a criterion for entanglement between the modes 
a+i and a_i (a state with this property is also referred 
as a two-mode squeezed state, or a continuous variable 
entangled state [8,20]). For later applications in the spin- 
1 condensate, we are more interested in two effective 

bosonic modes a'j_i, defined as a'^^i — a±iao/i/(aJao 

The modes a'^.^ are approximately bosonic if nearly all 
the atomic populations remain in the component oq, i.e., 
( ajflo ) ~ iV. The definition of the new modes a'^i, which 
are proportional to aj-i, is to remove the phase random- 
ness in the condensate by phase locking both of the modes 
a±i to the condensate mode oq. For the modes a'j.i, 
the quadrature phase amplitudes defined above reduce 
to /{AXiy\ ~ e±/2. The inequality ^l + f_ < 2 from 
the spin squeezing parameters thus becomes a sufficient 
criterion for entanglement between the modes a'_j_]^. As a 
result, a measurement of the spin squeezing parameters 
^^ provides criteria for both many-particle entanglement 
and two- mode entanglement. 

We will see that three-mode entanglement is also pro- 
duced in a spin-1 condensate. A general criterion for 
the three-mode entanglement is still absent, but our sys- 
tem to a good approximation will be in a pure three- 
mode state \4>)ii _i Q, and in this case the three- mode 
entanglement can be easily characterized by the quan- 
tity E3 ^ {/E{p+i)E{p^i)E{po), where E (p^) = 
—tr (pa log2 pa) denotes the von Neuman entropy of pa, 
and Pa is the reduced density operator of the mode a 
(a — +1,-1,0). It is obvious that for pure states E3. is 
non-increasing under local operations and positive if and 
only if the system state is intrinsically three-mode en- 
tangled. It will be shown that in a spin-1 condensate £"3 
can be significantly larger than 0. Note that different to 
the spin squeezing parameters it is not easy to directly 
measure Ej,. 

Now we turn to the physical system, which is a dilute 
gas of trapped bosonic atoms with hyperfine spin / = 1 . 
Assume that the system is very cold so that the collision 
interaction between atoms is effectively described by a 
pair-wise pseudo (5-potential, which preserves the hyper- 
fine spin of individual atoms and is rotationally invariant 
in the hyperfine spin space [12]. With this symmetry as- 
sumption, the most general form of the Hamiltonian can 
be dived into two parts H = Hg + Ha, with the symmet- 
ric and asymmetric Hamiltonians Hg and Ha respectively 
given by [12,15] (setting h = 1) 
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definition we have L^ = 2JjJ: and L^ = 2J^. Therefore, 
the Hamiltonian Ha expressed by J± has the form 



where ^a {a ~ — 1,0,+1) is the atomic field operator as- 
sociated with atoms in the spin state \f = i,mf — a). 
The first and the second terms of the Hamiltonian Hs 
denote the kinetic energy and the potential energy, with 
TO, the mass of the atom, and V, the trap potential, 
which has been assumed to be the same for all the 
three components. The 3x3 spin matrices S in Ha de- 
note the conventional 3-dimensional representation (cor- 
responding to the spin / = 1) of the angular momen- 
tum operator, with S*;^^ = {Sa,f3-i + Sa.,(3+i) /V^, S^^ = 
i {6a,f3-i — da,fi+i) /\/2, S^0 = aSaf3- The interaction en- 
ergy has been divided into the symmetric and the asym- 
metric parts, represented respectively by the last term of 
Hs and the Hamiltonian Ha- The symmetric interaction 
remains the same if we exchange two arbitrary spin com- 
ponents, whereas the asymmetric one is spin-dependent. 
The symmetric and the asymmetric collision coefficients 
As and Xa are determined respectively by the average 
and the difference of the scattering lengths in two differ- 
ent collision channels. For sodium atoms, Xa is positive 
and about 25 times smaller than A,, [11,12]. Thus, the 
Hamiltonian Hs is the dominant one for determining the 
dynamics of the motional state. However, the internal 
spin dynamics which we are interested in is completely 
determined by Ha since the symmetric part Hg is spin- 
independent. 

For generating entanglement, we start with a BEC 
prepared in the internal level \mf — 0). To quantify 
the amount of various kinds of squeezing and entangle- 
ment obtainable in this system, we need to calculate the 
time evolution of the entanglement criteria ^^ and E^ 
defined before. For this purpose, we make use of the 
single motional mode approximation to the Hamiltonian 
Ha following Refs. [15-17]. Note that the Hamiltonian 
Hs is the dominant one for determining the motional 
state. When the system is at a very low temperature, 
the motional state is frozen to be approximately the 
ground state ip (r) (normalized a,s J dr \ip {r)\ = 1) of the 
spin-independent Hamiltonian Hs, and we can factorize 
the field operators as ^q, « aa(p{r) {a = 0,±1), where 
tta is the usual bosonic annihilation operator. With 
this approximation, one can introduce a spin operator L 
with L = J2a B (^a^ a 13 a (3, and the spin-dependent Hamil- 
tonian Ha is simplified to Ha = X'a (L^ — 2N) , where 
A^ = Xa J dr \(f (r)| /2. This Hamiltonian has been used 
for studying the spin ground state structure [15]. We 
use it here for calculating squeezing and entanglement in 
this system. For this purpose, it is convenient to express 
the Hamiltonian Ha by the more relevant spin operators 
3± defined before. For the specific system configuration 
considered here, L^ — for the initial state, and it will 
remain zero since L^ is conserved. On the other hand, by 



Ha = 4A1 (jf + Jl^ 
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where we have dropped the irrelevant constant 2A^iV. 
The only non-zero commutation relation between J+ and 

J- is [JlJl] = [Jl,Jl\ - (aLia+i-a^ia-i)/4, 
which remains very small if the population of the atoms 
is still dominantly in the level \mf = 0), as it is in our 
scheme. In this case, as an approximation we can con- 
sider separately the dynamics of J+ and J_, respectively 
with the effective Hamiltonian J^^ or J^ . It is well 
known that the Hamiltonian J^ or J^ will produce a 
spin squeezing state, with the minimum squeezing pa- 
rameter about (3/iV) /2 after an evolution time pro- 
portional to N~'^/^ /X'a [18]. This is the intuitive rea- 
son why we can get substantial spin squeezing and en- 
tanglement in this system. To be more quantitative, 
we have calculated exactly the time evolution of the 
squeezing parameter ^^ under the Hamiltonian (5) us- 
ing the numerical approach. The system spin state \4>) 
under the Hamiltonian (5) can be expanded as \4>) = 

X^n *-" 1-^ ~ 2'^' "^j '^)o +1 -1 ^^ ^^^ number basis, where 
the subscript a [a — Q, ±1) represents the corresponding 
internal mode a^. One can solve numerically the time 
evolution of the expansion coefficients c„, from which the 
spin squeezing parameters ^^ are exactly obtainable. We 
find that the squeezing magnitudes ^± = ming ^^ evolve 
in the same way for 3^ and J , while their squeezing di- 
rections corresponding to the optimal angle 6 are always 
orthogonal to each other. This means that the criterion 
mine i£,+ +^-^ ) /2 for entanglement between the two 

effective modes a'j_i evolves in the same way as the crite- 
rion ^± for many-particle entanglement . The resultant 
time evolution of ^± is shown in Fig. la, together with 
a prediction from the simple Hamiltonian J^^ (or J^ ). 
One can see that the results from the two approaches 
agree quite well, with only slight difference in the maxi- 
mum obtainable squeezing and the corresponding evolu- 
tion time. For 10^ atoms, about three orders reduction is 
obtainable in S,± , which is a clear demonstration that sub- 
stantial many-particle entanglement and two-mode en- 
tanglement have been generated. Fig. lb shows the time 
evolution of the three- mode entanglement E^. One can 
see that significant three-mode entanglement can be gen- 
erated with only 10"^ atoms. 

Finally, we would like to address some imperfections 
and approximations used in this scheme. First, in a prac- 
tical case there would be some remaining external mag- 
netic field due to the special way of preparing the con- 
densate. It is remarkable that the linear Zeeman shift 
from a known homogeneous magnetic field has no in- 
fluence on our scheme since it commutes with the spin- 
dependent Hamiltonian Ha and will be eliminated if one 



transforms to the interaction picture. As a result, only 
quadratic Zeeman shift and linear Zeeman shift from an 
inhomogeneous field have influence. These high-order ef- 
fects can be easily controlled very small in the relevant 
time scale. For instance, one can see from Fig. 1 that 
the time scale in this scheme is roughly set by the inverse 
of iVAJj, which is about 27r x lO^Hz for a typical density 
3 X lO^'^cm"'^ of sodium atoms [11]. Within this time, 
the quadratic and the inhomogeneous Zeeman shifts have 
negligible effects if the magnetic field and its gradient 
are respectively smaller than ImG and ImG/cm. Sec- 
ondly, we have assumed that the three components with 
ruf = 0, ±1 have the same spatial wave function. This 
is a good approximation as long as there is no demixing 
instability leading to phase separation of different com- 
ponents [2,22]. For the two-component condensate with 
components a and /3 having equal macroscopic popula- 
tions, a mean-field theory has shown that there is no 
demixing instability if the intcr-componcnt collision co- 
efficient Xa/3 is smaller than the geometric average of the 
intra-component collision coefficients ^JXaaXpii [22,11]. 
For the relevant dynamics in our scheme with only the 
level |m/ = 0) macroscopically populated, Fig. 1 shows 
that this system is well approximated by two independent 
sets of two-component condensates. For each set of two- 
component condensate, we can apply the above result, 
and find there is no demixing instability in this scheme 
if the asymmetric collision coefficient Aq is positive, as it 
is for sodium atoms. This is different to the case where 
both of the levels with mj — 0, -1-1 are equally populated, 
which leads to demixing instabilities [10,11]. Finally, we 
assumed the single-mode approximation for the atomic 
motional state. This approximation is supported by the 
following arguments: first, in the two-component conden- 
sate more complete numerical simulations gives a squeez- 
ing which agrees quite well with the calculation from the 
single-mode approximation [2,6,23], and the present spin- 
1 system is well approximated by two independent sets of 
two-component condensates; and second, a direct energy 
analyses in Ref. [16] also shows the validity of the single- 
mode approximation for the spin-1 condensate with only 
the level |m/ = 0) initially populated. 

In summary, we have shown how to directly generate 
and detect various kinds of squeezing and entanglement 
in the experimentally realized spin-1 condensate. This 
helps for the final observation of quantum entanglement 
in Bose-Einstein condensates. 

This work was supported by the Austrian Science 
Foundation, the EQUIP, the ESF, the European TMR 
network Quantum Information, and the Institute for 
Quantum Information Gcs.m.b.H. L.M.D. thanks in ad- 
dition the support from the Chinese Science Foundation. 



[i; 



[7] 
[8] 

[9; 

[lo; 
[11 



M. H. Anderson et al, Science 269, 198 (1995); Davis, 
K.B. et al, Phys. Rev. Lett. 75, 3969 (1995); C. C. 
Bradley et al, Phys.Rev. lett. 75, 1687 (1995). 
A. S0rensen, L.-M. Duan, J. I. Cirac, and P. Zoller, Na- 
ture 409, 63 (2001). 

H. Pu and P. Meystre, Phys. Rev. Lett. 85, 3987 (2000). 
L.-M. Duan, A. S0rensen, J. I. Cirac, and P. Zoller, Phys. 
Rev. Lett. 85, 3991 (2000). 
N. P. Bigelow, Nature 409, 27 (2001) 



U. Poulsen a nd K. M0lmer, ^uant-ph/0101089 



mat/0006030 



cond- 



K. Helmerson and L. You, quant-ph/010413C 



Special issue on quantum information. Phys. World 11, 

33 (1998). 

D. M. Stamper-Kurn et al, Phys. Rev. Lett. 80, 2027 

(1998). 

J. Stenger et a/.. Nature 396, 345 (1998). 

D. M. Stam per-Kurn and W. Ketterle, 



cond-mat/0005001 



[12] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[13] T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 
(1998). 

[14] E. V. Goldstein and P. Meystre, Phys. Rev. A 59, 1509 
(1999). 

[15] C. K. Law, H. Pu, and N. P. Bigelow, Phys. Rev. Lett. 
81, 5257 (1998). 

[16] H. Pu et al, Phys. Rev. A 60, 1463 (1999). 

[17] T.-L. Ho and S. K. Yip, Phys. Rev. Lett. 84, 4031 (2000). 

[18] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993). 

[19] D. Wineland, J. J. Bollinger, W. M. Itano, and D. J. 
Heinzen, Phys. Rev. A 50, 67 (1994). 

[20] L.-M. Duan, G. Giedke, J. I. Chac, and P. Zoller, Phys. 
Rev. Lett. 84, 2722 (2000). 

[21] For instance, in the second 

quantization, the state |A/2, A/2) , ^ _-^ is particle en- 
tangled, but not mode entangled; and the single-particle 
state (|1,0) I j^ _^ -I- |1,0) I j^ _jj /-\/2 is mode entangled, 
but not particle entangled. 

[22] A. Sinatra and Y. Castin, Eur. Phys. J. D 8, 319 (2000). 

[23] D. Jaksch, private communication. 




la lb 

FIG. 1. a. Time evolution of the squeezing magnitudes for 
10^ atoms. The solid curve shows the squeezing from an exact 
numerical solution to the Hamiltonian (5), and the dashed 
curve is the result from the Hamiltonian 4AJ, J^^. The upper 
dotted curve shows the ratio of the atomic population in the 
level |m/ = 0), which is basically 1. b. Time evolution of the 
three-mode entanglement £3 for 10"^ atoms. 



